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ABSTRACT 

-xj- : 

s ! ■ I discuss the impact of color-octet contributions to inclusive J/tp photoproduc- 

tion at HERA. Emphasis is put on resolved photon processes and on J/t/j po- 
larization, which will be experimentally accessible at HERA in the near future. 

Qs, | Both analyses provide a powerful test of the NRQCD factorization approach to 

charmonium production. 

-xf ' 

^ 1. Introduction 

t> 

The production of heavy quarkonium states in high-energy collisions provides an 
important tool to study the interplay between perturbative and non-perturbative 
QCD dynamics. While the creation of heavy quarks in a hard scattering process can 
be calculated in perturbative QCD, the subsequent transition to a physical bound 
state introduces non-perturbative aspects. A rigorous framework for treating quarko- 
nium production and decays has been developed only recently.0 The factorization 



approach is based on the use of non-relativistic QCD (NRQCD )□ to separate the 
short- distance parts from the long-distance matrix elements and explicitly takes into 
account the complete structure of the quarkonium Fock space. According to this 
formalism, the inclusive cross section for the production of a quarkonium state H can 
be expressed as a sum of terms, each of which factors into a short-distance coefficient 
and a long-distance matrix element: 

da(H + X) = J2MQQ[n}+X)(0 H [n}} (1) 

n 

Here, da denotes the short- distance cross section for producing an on-shell QQ-p&ii 
in a colour, spin and angular-momentum state labelled by n. The universal NRQCD 
matrix elements (O h [n]) give the probability for a QQ-p&ir in the state n to form 
the quarkonium state H . The relative importance of the various terms in ([]]) can be 
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estimated by using NRQCD velocity scaling rules.0 For v — ► (v being the average 
velocity of the heavy quark in the quarkonium rest frame) each of the NRQCD matrix 
elements scales with a definite power of v and the general expression (JI]) can be 
organized into an expansion in powers of the heavy quark velocity. 

The general factorization formula (|I[) is believed to hold also for quarkonium pro- 
duction in hadron-hadron or photon-hadron collisions. In the case of inclusive open 
heavy flavour production it has been shown that the hadronic cross section factorizes 
into a partonic hard scattering cross section multiplied by light quark and gluon par- 
ton densities, with higher-twist corrections being suppressed by powers Aq CD /raQ.S 
The argument involves averaging over a sufficiently large range of the heavy quark 
transverse momentum and uses the fact that the main contribution to the total 
inclusive cross section comes from p T ~ rriq ^> Aq CD . Unfortunately, the situation 
is more involved for quarkonia production in hadronic collisions, since the region of 
small transverse QQ momentum significantly contributes to the total cross section. At 
small transverse momentum however, the intermediate colour-singlet or colour-octet 
QQ pair moves parallel with a remnant jet and might interact before the physical 
bound state has been formed. Only in the heavy quark limit, where all scales in- 
volved in the bound state physics are much larger than Aqcd, higher- twist terms are 
necessarily strongly suppressed. For charmonium production in hadronic collisions, 
it seems however reasonable to restrict the analysis to the region px 3> Aqcd, where 
the general factorization formula should safely be applicable. Excluding the small 
Pt region is also necessary to guarantee perturbative stability in higher-order QCD 
calculations of J/ip photoproduction.ifl 

The NRQCD formalism implies that colour-octet processes, in which the heavy- 
quark antiquark pair is produced at short distances in a colour-octet state and sub- 
sequently evolves nonperturbatively into a physical quarkonium, must contribute to 
the cross section. As discussed extensively in the literaturei§, colour-octet contribu- 
tions appear as the most plausible explanation of the large direct ip production cross 
section observed at the Tevatron0. The NRQCD approach is certainly the correct 
theory for quarkonium production in the heavy quark limit. It is however not clear 
whether the charm quark mass is sufficiently large to allow for a reliable expansion in 
the heavy quark velocity. In order to establish the phenomenological significance of 
the colour-octet mechanism and the universality of the NRQCD matrix elements for 
charmonium production it is thus necessary to identify colour-octet contributions in 
different production processes. The analysis of J/ip photoproduction at HERA ap- 
pears to be a very powerful tool to constrain the colour-octet matrix elements and to 
test the picture of quarkonium production as developed in the context of the NRQCD 
factorization approach. 

The production of J/ip particles in high energy ep collisions at HERA is dom- 
inated by photoproduction events where the electron is scattered by a small angle 

"Higher- twist corrections to colour-singlet J/ ip photoproduction have been analysed recently.! 



producing photons of almost zero virtuality. The measurements at HERA provide 
information on the dynamics of J/ip photoproduction in a wide kinematical region, 
30 GeV y/s 1P <j 200 GeV, corresponding to initial photon energies in a fixed-target 
experiment of 450 GeV ^ E 1 ^ 20,000 GeV. The production of J/ip particles in 
photon-proton collisions proceeds predominantly through photon-gluon fusion, where 
the photon interacts directly with the gluon from the proton. Besides the direct 
photoproduction channel, J / ip production at HERA can also take place via resolved 
photon contributions where the photon behaves as a source of partons which interact 
with the partons in the proton. Resolved processes are expected to contribute sig- 
nificantly to the lower endpoint of the J / ip energy spectrum and might be probed at 
HERA for the first time in the near future. 6 



2. Direct Photon Contributions 



For J/ip production and at leading order in the velocity expansion, the general 
factorization formula ([[]) reduces to the standard expression of the colour-singlet 
mo delH. The short-distance cross section for J/ip photoproduction through direct 
photons is given by the subprocess 

l + g^cc[ 3 S u l}+g (2) 

with cc in a colour- singlet state (denoted by I), with zero relative velocity, and 
spin/angular-momentum quantum numbers 2S+1 Lj = 3 S\. Relativistic corrections 
to the colour-singlet channel due to the motion of the charm quarks in the J/ip 
bound state enhance the small-p t region, but can be neglected for p^ <; 1 GeV. Ill The 
calculation of the higher-order perturbative QCD corrections to the short-distance 
cross section (0) has been performed recent lyffl Inclusion of the NLO corrections 
reduces the scale dependence of the theoretical prediction and increases the cross sec- 
tion significantly, depending in detail on the 7p energy and the choice of parameters.0 
Details of the calculation and a comprehensive analysis of total cross sections and 
differential distributions for the energy range of the fixed-target experiments and for 
J/ip photoproduction at HERA can be found elsewhere.i 

The leading colour-octet configurations which contribute to J / ip photoproduction 
at pt 7^ are produced through the subprocesses0 

7 + 5 - cc[ 1 ^o, 3 5i, 3 P ,i,2,8] + ( ? 
1 + q(q) - cc[ 1 5 , 3 5 1 , 3 P ,i,2,8] + g(g) (3) 

The transition of the colour-octet cc [ 2S+1 Lj, 8] pairs into a physical J/ip state through 
the emission of non-perturbative gluons is described by the long-distance matrix 

^Elastic mechanisms only contribute to ±he region of small J/ip transverse momentum and can be 
suppressed by an appropriate cut in pr .til 



elements (O j I^ [ 2S+1 Lj, 8]). They have to be obtained from lattice simulations or 
measured directly in some production process. According to the velocity scaling 
rules of NRQCD, the colour-octet matrix elements associated with S-w&ve quarkonia 
should be suppressed by a factor of v 4 ~ 1CT 2 compared to the leading colour-singlet 
matrix element. However, explicit calculation of the processes (|3|) shows that the 
short- distance factors of the [3So,8] and [ 3 Po,2,8] channels are strongly enhanced as 
compared to the colour-singlet term and compensate the O(l0~ 2 ) suppression of the 
corresponding non-perturbative matrix elements. This short-distance enhancement is 
due to t-channel gluon exchange which contributes to the pSo, 3 Po,2, 8] cross sections, 
but not to the leading-order colour-singlet channel. 

3. Resolved photon contributions 

Inclusive J/if> photoproduction can also take place via resolved photon interac- 
tions, where the photon first splits into a flux of light quarks and gluons which then 
may fuse with a gluon or light quark from the proton to form a cc pair. This pro- 
cess has been extensively analyzed in the past in the framework of the colour-singlet 
model, where the cross section is given by gluon-gluon fusion into a 3 S± colour-singlet 
cc pair, which subsequently hadronizes into a J/if>: 

g + g^cci%,l}+g (4) 

The process contributes to the overall cross section only marginally except near 
the lower endpoint of the J /if) energy spectrum.lll 

Within the NRQCD factorization approach, however, more resolved-photon chan- 
nels have to be considered, where J/ if) production proceeds via colour-octet cc pairs.0 
The leading colour-octet contributions are: 

g + g - cc[ 1 5 , 3 S 1 , 3 P 0il , 2 ,8]+^ 
g + q(q) - cc [%, 3 S X , 3 P 0)1 , 2 , 8] + q(q) 

q + q -> cc[ 1 5o, 3 ^ 1 , 3 P ,i,2,8] + (7 (5) 

They are in every respect analogous to the ones which have been argued to strongly 
increase the J/if> production cross section in pp collisions at the Tevatron.i'i Colour- 
octet contributions to resolved J/ if> photoproduction thus provide an important direct 
test of the NRQCD explanation of the Tevatron data. 

4. The J /if) Energy Distribution 

The J /if} energy distribution da/dz offers unique possibilities to study the rela- 
tive importance of the different mechanisms that contribute to J/ if) photoproduction. 
The scaling variable z is defined by z — p • k^p / p • fc-y, p, k^,^ being the momenta of 



the proton and J/ip, 7 particles, respectively. In the proton rest frame, z is the ra- 
tio of J/tp to 7 energy, z = E^/E 7 . Adopting NRQCD matrix elements consistent 
with those extracted from the fits to prompt J/ip data at the Tevatroni (see Table 
[]]), one finds significant colour-octet contributions near the upper and lower end- 
point of the J/ip energy spectrum. This is demonstrated in Fig. [I], where I plot the 

Table 1. Values of the NRQCD matrix elements used in the numerical analysis, with the 
velocity and mass scaling. 
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colour-singlet and colour-octet contributions to da/dz at a typical HERA energy of 
y/s-yp = 100 GeV in the restricted range p± > 1 GeV via direct and resolved photon 
pr o cesses. c The colour-octet processes are predicted to exceed the colour-singlet ones 
by more than an order of magnitude for z J> 0.6 and z <, 0.2. Resolved contributions 
are negligible for z J> 0.3, but colour-octet terms are expected to become significant 
at the lower end of the energy spectrum. d The theoretical predictions are compared 
to experimental results obtained by the Hl@ and ZEUsH collaborations at HERA. 
From Fig. p] one can conclude that the large colour-octet enhancement in the region 
z ^ 0.6 is not supported by the experimental data and that the J/ip energy spectrum 
is adequately described by the colour-singlet channel. Taking into account the large 
uncertainty of the cross section normalization due to the choice for the charm quark 
mass and the QCD coupling, the current experimental data on the J/ip energy distri- 
bution can be accounted for by the colour- singlet channel alone. This can be judged 
by comparing the experimental data with the shaded band in Fig. [I] which indicates 
how much the colour- singlet cross section is altered if m c and A^ 4 -* vary in the range 
1.35 GeV <m c < 1.55 GeV and 200 MeV < A< 4 > < 300 MeV. e 

One should be careful to interpret the discrepancy between the theoretical pre- 
dictions for the colour-octet contributions and the HERA data shown in Fig. [I] as a 
failure of the factorization approach. As pointed out recently, the NRQCD expansion 
breaks down close to the phase space boundary z — > 1, and no prediction can be 
made for z <; 0.75 without resumming large corrections in the velocity expansion.!^ 
Therefore, the shape of da/dz near the endpoint region cannot be used to access 
the importance of colour-octet contributions. Nevertheless, the analysis of the the 

c Note that a cut in the J/tp transverse momentum, pt > 1 GeV, is sufficient to exclude the elastic 
contributions; no additional cut in z is required. 

"Similarly, colour-octet processes enhance the \ow-z region in J/ip photoproduction at large px via 
fragmentation mechanisms .EJ !_. 
e To study the a s dependence of the cross section, I use consistently adjusted sets of parton densities.^] 
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Fig. 1. Colour-singlet (CS) and colour-octet (CO) contributions to the J/ip energy distribution 
da/dz at the photon-proton centre of .mass energy s lp = 100 GeV integrated in the range pr > 
1 GeV, compared to experimental dataliao. Both direct and resolved photon processes are included. 
Parameters: m c = 1.5 GeV, renormalization/factorization scale fi — 2m c , GRV parton distribution 
functionsEJ with A' 4 ) = 200 MeV. The shaded band reflects the uncertainty in the normalization of 
the direct colour-singlet contributions due to variation of m c and a s , as described in the text. 

J/ip energy distribution provides a crucial constraint on the colour-octet matrix el- 
ements if one averages over a sufficiently large interval Az 3> v 2 ~ 0.25 containing 
the endpoint z = 1, as shown in Fig. |2|. Colour-octet mechanisms should contribute 
significantly to the region 0.6 ^ z < 1 and exceed the colour-singlet cross section by 
more than an order of magnitude. While the overall size of the cross section is subject 
to large uncertainties due to m c and a s , the distinctive modification of the J/ip energy 
spectrum by colour-octet processes in the low and large z region should clearly be 
visible in the experimental data. Moreover, the low-z region of the energy spectrum 
is not expected to be sensitive to higher order terms in the velocity expansion and 
can thus reliably be predicted within leading-twist NRQCD. 

The results at HERA seem to indicate that the values of the colour-octet matrix 
elements (O j ^ [3So,8]) and (O j ^ [ 3 Po,S\) are considerably smaller than suggested by 
the fits to the CDF data at the Tevatron. This does however not necessarily imply the 
non- universality of the NRQCD matrix elements for J/ip production. The determina- 
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Fig. 2. Same as Fig. |l], averaged over large intervals in z. 

tion of the (O j ^ [ 1 So> 8]) and (O j I^ [ 3 Pq, 8]) matrix elements from the Tevatron data 
is very sensitive to all effects that modify the shape of the J / ip transverse momentum 
distribution at pt ^ 5 GeV and might yield much too large values.il The theoretical 
uncertainties include the small-x behaviour of the gluon distribution, the evolution of 
the strong coupling, next-to-leading order contributions, intrinsic transverse momen- 
tum of the partons in the proton, and also systematic effects inherent in NRQCD, 
such as the inaccurate treatment of the energy conservation in the hadronization of 
the colour-octet cc pairs. With higher statistics data at HERA it will be possible to 
extract more detailed information on the colour-octet matrix elements, in particular 
from the analysis of the J/if> energy distribution. 



5. J/ip Polarization 

The NRQCD factorization approach allows for unambiguous predictions of the 
quarkonium polarization.0 Transverse polarization of J / if) and ip' particles, produced 
directly at large transverse momentum in pp collisions at the Fermilab Tevatron, has 
in fact emerged as the most prominent test of colour-octet contributions and spin sym- 
metry in charmonium production.^ Polarization signatures in J/ip photoproduction 



have been analysed in detail in the past within the colour- singlet model.c3 Recently, 
a comprehensive calculation has been performed in the context of the NRQCD fac- 
torization approach, including colour- singlet and colour-octet contributions to both 
direct as well as resolved photon processes.^ Results are available for the most gen- 
eral decay angular distribution in any given reference frame. Here, I only consider the 
polar-angle distribution in the s-channel helicity frame, in which the polarization axis 
in the J/ip rest frame is defined as the direction of the J/ip three-momentum in the 
photon-proton cms frame. The polar angular distribution in the decay J/ip — > 
is given by 

oc 1 + a cos 2 6, (6) 

dcosfl 

with 9 the angle between the lepton three- momentum in the J/ip rest frame and the 
polarization axis. Fig. |3| displays colour- singlet and colour-octet contributions to the 
polar asymmetry a, defined in (^), as function of the J/ip energy variable z. It is worth 
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Fig. 3. Colour-singlet (CS) and colour-octet (CO) contributions to the polar angle asymmetry 
a in (^) at the photon-proton centre of mass energy ^s lp = 100 GeV integrated in the range 
pr > 1 GeV. Both direct and resolved photon processes are included. 

pointing out that the asymmetry predictions involve ratios of cross sections and are 
thus to a large extent insensitive to uncertainties in the charm quark mass, the QCD 
coupling or parton distribution functions. Therefore, polarization measurements are a 



rather clean test of the underlying quarkonium production mechanisms and will allow 
to access the relative importance of colour-singlet and colour-octet contributions. 



6. Conclusions 

I have discussed colour-singlet and colour-octet contributions to J/ip photopro- 
duction at HERA. The high-statistics data to be expected in the near future will 
allow the analysis of resolved photon processes and J/ ip polarization to further test 
the phenomenological significance of colour-octet contributions. Various other chan- 
nels and final states can be accessed in the future at HERA, likeohotoproduction 
of ip', T and Xc particleJiE, associated J/ip + 7 productionlil'il'ii, fragmentation 
contributions at large as well as deep inelastic J /if) production^. The future 

analyses of inclusive quarkonium production at HERA offer unique possibilities to 
test the NRQCD factorization approach to charmonium production. 
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